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Abstract 

This is a paper on young jet material in a frustratingly complex environment. 

NGC 1052 has a compact, flat- or GHz-peaked-spectrum radio nucleus con- 
sisting of bi-symmetric jets, oriented close to the plane of the sky. Many features 
on both sides move away at u app ~ 0.26c (Ho = 65kms~ 1 Mpc -1 ). VLBI at 
seven frequencies shows a wide range of spectral shapes and brightness temper- 
atures; there is clearly free-free absorption, probably together with synchrotron 
self-absorption, on both sides of the core. The absorbing structure is likely to be 
geometrically thick and oriented roughly orthogonal to the jets, but it is patchy. 

Hi VLBI shows atomic gas in front of the approaching as well as the receding 
jet. There appear to be three velocity systems, at least two of which are local to 
the AGN environment. The "high velocity system", 125 to 200 km s -1 redward 
of systemic, seems restricted to a shell 1 to 2 pc away from the core. Closer to the 
centre, this gas might be largely ionised; it could cause the free-free absorption. 

WSRT spectroscopy shows 1667 and 1665 MHz OH absorption over a wide 
velocity range. OH and Hi profile similarity suggests co-location of molecular and 
atomic "high velocity" gas; the connection to H2O masing gas is unclear. Further, 
at "high velocity" we detected the OH 1612 MHz satellite line in absorption and 
the 1720 MHz line in emission, with complementary strengths. 
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1 Introduction 



The nearby 1 LINER (e.g., Gabel et al. 2000) elliptical galaxy NGC 1052 has an unusu- 
ally prominent central radio source (1 to 2 Jy). It is variable on timescales of months 
to years (e.g. Heeschen & Puschell 1983) and has a fairly flat spectrum between 1 and 
30 GHz, which has sometimes been classified as Gigahertz peaked (e.g., de Vries et al. 
1997). The overall radio structure is core-dominated, and has two lobes spanning only 
about 3kpc (Wrobel 1984), so that NGC 1052 meets the traditional size limit for CSSs, 
but not for CSOs. With VLBI, detailed sub-parsec scale scrutiny of the active nucleus 
and its inner environment are possible. 

A more extensive analysis of the data presented below is in Vermeulen et al. (2002). 

2 Kinematics 

Ten epochs of 15 GHz VLBA data (Fig. 1) show a two-sided source, with oppositely 
directed, slightly curved jets and a prominent gap 0.1 to 0.2 pc west of the brightest fea- 
ture in most images. At our high linear resolution, NGC 1052 shows complex evolution, 
which we have quantified using a consistent set of moving components, rank-numbered 
"E"ast and "West by decreasing distance to the centre (Fig. 1). While not every one 
of these should be interpreted as a true physical entity (a plasmon or a "cannon-ball"), 
we are confident that the set as a whole gives an appropriate representation of the 
motions in the jets. There is a plausible relative alignment of the epochs, in which the 
motions of most of the components, on both sides of the centre, are consistent with 
being ballistic: linear on the sky and constant over time. Features on the two sides 
move in opposite directions with roughly equal apparent velocities of 0.26 ± 0.04 c. 
Using these error margins, the jets are oriented at most 33° from the plane of the sky. 

3 Ionised gas 

The left panel in Fig. 2 shows data at 43, 22, 15, 8, 5, 1.6, and 1.4 GHz from VLBA 
observations within a few days of each other in July 1997. There is a wide range of 
spectral shapes, which proceed from steep, through convex, to highly inverted, from 
the outer jets towards the middle, and produce a distinctive central hole in the VLBA 
images. On both sides of the gap there are components with a low-frequency brightness 
temperature well below 10 10 K, and a low-frequency spectral cutoff steeper than a = 3 
if expressed as a power law. As first reported in Kellermann et al. (1999), we find that 
the only plausible explanation is free-free absorption from ionised gas; a conclusion also 
reached with other data by Kameno et al. (2001) and Kadler et al. (2001a). At about 
0.5 pc along the eastern jet, and 1 to 2pc along the western jet, the signature of free- 
free absorption becomes hard to distinguish from synchrotron-self-absorbed spectral 

1 We adopt Ho=65kms _1 Mpc -1 , with no corrections for local deviations from the Hubble flow, so 
that the optical stellar absorption line heliocentric redshift, cz = 1474kms~ 1 (Sargent et al. 1977), 
corresponds to a distance of 22 Mpc and implies that 1 mas « 0.1 pc 
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Figure 1: Left: VLBA 15 GHz contour images of NGC1052 at ten epochs, aligned 
following the kinematic analysis of Sect. 2, and shown spaced by their relative time 
intervals. All contour levels increase by factors of 3 1 / 2 , starting at 0.58% of the peak 
brightness in each image. The restoring beams are all 0.1x0. 05 pc in P.A. 0°. 
Right: Distances to the adopted origin for the Gaussian components (see Sect. 2) fitted 
to these datasets. Lines show the best-fit linear velocities. 



shapes. If conditions are such that similar peak frequencies result, the two absorption 
mechanisms are difficult to disentangle, but it is likely that both play a role. 

The absorption is more pronounced along the western jet than at corresponding 
locations along the eastern jet. This suggests a geometrically thick disk- or torus-like 
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Figure 2: Left: All flux densities from a multi-frequency multi-component VLBI anal- 
ysis. The diversity of spectral shapes is emphasised by the frequency colour coding. 
Right: Flux density changes with time for some Gaussian components fitted to the 
multi-epoch 15 GHz observations shown as a function of their distance to the origin. 



absorbing region, more or less perpendicular to the jets, which are oriented close to 
the plane of the sky, with the eastern jet approaching and the western jet receding. 
Variations on this scenario with a thinner disk or torus, either not oriented orthogonal 
to the jets, or warped, are possible. Furthermore, the range of spectral shapes and 
the flux density variations of components tracked over time at 15 GHz (Fig. 2, right 
panel), suggest that the absorbing region has a fairly well-defined overall geometry but 
substantial patchiness in detail. It is then not very meaningful to use a single epoch 
for detailed modeling. The deepest absorption seen, r ~ 1 at 43 GHz over the central 
region, implies a volume density of n e ~ 10 5 cm -3 if the free-free absorbing gas were 
distributed uniformly along a path-length of 0.5 pc with a temperature T = 10 4 K. 



4 Atomic gas 

Figure 3 depicts Hi VLBI spectra from July 1998 at all locations where the signal- 
to- noise ratio allowed detection of absorption at the level of a few percent. We think 
there are three different absorption systems, at least two of which are probably due to 
atomic gas on parsec or sub-parsec scales, local to the AGN environment. 

The most remarkable system is at "high velocity", redshifted by 125 to 200kms _1 
with respect to the stellar systemic velocity. It gives rise to 5 to 20 % absorption at 1 
to 2pc along both the western, receding jet, and the eastern, approaching jet, with a 
possible west-to-east velocity gradient of some 10kms~ 1 pc~ 1 . But the high velocity 
system is absent closer to the core, at least on the eastern side, where it could easily 
have been detected. The combination of free-free absorption covering the inner parsec, 
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Figure 3: Left: Hi optical depths observed in 1998 at various locations along the jets 
of NGC 1052. Multiple offsets are used for clarity; dotted lines show the zero levels. 
Right:WSRT optical depth measurements (negative is emission) for all four 18 cm OH 
lines. The offset zero lines are indicated; also note the different scale for our integrated 
VLBI Hi spectrum, shown convolved to the same resolution for comparison. 



and atomic gas in an annulus at 1 to 2 pc, is quite natural: the innermost region and/or 
the surface of an accretion disk or torus receive the most intense ionising radiation. 
Weaver et al. (1999), Guainazzi & Antonelli (1999), Guainazzi et al. (2000), and Kadler 
et al. (2002b) discuss soft X-ray absorbing gas towards the nuclear continuum source, 
with large column depths (iV H = 10 23 to 10 24 cm -2 ), and possibly a patchy distribution, 
or involving two components with a substantially different density, which matches the 
evidence we see for patchiness in the radio spectra and flux density evolution of jet 
components. An Hi optical depth of 20% with a FWHM of 20kms -1 implies a column 
depth of Nn = 10 21 T SPi ioo cm -2 , but close to an AGN the spin temperature T sp may 
well be one or two orders of magnitude above 100 K (e.g., Maloney et al. 1996). 

5 Molecular gas 

Figure 3 shows WSRT spectra of the 18 cm OH lines. We find that 1667 and 1665 MHz 
absorption, recently detected by Omar et al. (2002), extends over at least as wide a 
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velocity range as Hi. The peak 1667 MHz depth, 0.4% in the high velocity system, 
suggests a column depth of order 10 14 cm~ 2 . The 1667/1665 ratio ranges from near 1 
at low velocities to approximately 2 in the high velocity system. We have discoverd 
that the satellite lines are also present: 1612 MHz in absorption and 1720 MHz in 
emission. Their conjugate profiles probably result from excitation in a far infra-red 
radiation field when the OH column density is sufficiently large; competing pumping 
mechanisms determine which line is in emission and which one in absorption in specific 
density and temperature regimes, as modeled for Cen A by van Langevelde et al. (1995). 

The OH main lines and the total Hi profile are remarkably similar in the high ve- 
locity system, and we suggest co-location of these high velocity atomic and molecular 
gas components. They probably do not coincide with the H 2 masers at 0.1 to 0.2 pc 
along the receding jet (Claussen et al. 1998), even though these are at the same ve- 
locity. Other questions also remain regarding the nature of the high velocity system. 
Interpretation of the velocity gradient in Hi as evidence for a structure rotating around 
the nucleus is contradicted by the fact that the centroid is redshifted by 150 km s -1 or 
more from the systemic velocity. But if it is instead infalling gas, the nature of the 
central hole in Hi is unclear. 
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